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20.  Abstract  (cont'd) 


Experimental  data  verify  this  conclusion,  as  agreement  with  the  calculated 
value  of  power  iv  better  than  1%.  Losses  due  to  eddy  currents  and  hysteresis 
are  found  to  be  much  less  serious  amount ! ng  to  no  more  than  a few  percent, 
fteconmendat ions  of  possible  remedies  for  the  saturation  problem  are  made. 

A computer  study  is  made  to  determine  the  effect  of  variation  of  the  relevant 
design  parameters  on  the  peak  voltage  output  of  the  generator.  The  resulting 
design  matrices  provide  a useful  guide  to  design  optimization  as  well  as  a 
clear  delineation  between  the  parameter  combinations  which  result  in  saturation 
of  the  reed  material  and  those  which  do  not.  Aeconaendat ions  are  made  for 
design  changes  to  improve  the  generator  performance. 
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A FLUX  CIRCUIT  ANALYSIS  FOR  THE  MAGNETIC  TRANSDUCER 
OF  A FLUIDIC  REED  GENERATOR 


INTRODUCTION 

There  is  considerable  interest  in  the  utilization  of  fluidic  generators 
in  various  fuzing  applications.'*^  These  devices  are  mounted  in  the  ogives 
of  various  projectiles  and  they  essentially  convert  some  of  the  energy  of 
the  relative  motion  of  projectile  and  air  first  into  acoustical  energy  and 
thence  to  electrical  energy.  A typical  device  is  pictured  schematically  in 
Fig.  I. 
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Fig.  I,  Schematic  View  of  Fluidic  Generator  with  Reed  Type  Magnetic  Transducer 


I.  Carl  J.  Camtagnoulo,  ‘The  Fluidic  Generator,"  HOI  Technical  Report-1328, 
September  IJ66. 

i,  Carl  J,  Canpagnoulo,  "Fluidic  Power  Generators  for  Ordnance  Application," 
MDL  Technical  Report-1423*  December  1968. 

3.  Carl  J.  Campagnoulo  and  Richard  N.  Gottron,  "The  Fluidic  Generator:  A 
New  Electrical  Power  Source,"  24th  Annual  Proceedings  Power  Sources 
Symposium,  May  19-21,  1970. 
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As  the  projectile  moves  forward,  air  is  forced  through  the  annular  orifice 
to  excite  the  resonant  cavity  beyond.  The  oscillations  of  the  cavity  in 
turn  cause  the  diaphragm  to  vibrate.  This  vibration  is  then  transmitted  to 
a permeable  reed  of  Armco  electrical  steel.  As  can  be  seen  from  Fig.  I, 
when  the  reed  is  in  its  equilibrium  position  it  is  exactly  midway  between 
the  pole  pieces  of  both  magnetic  gaps  and,  therefore,  carries  no  magnetic 
flux.  The  vibration  of  the  reed,  however,  changes  the  relative  lengths  of 
gaps  I and  2 and  flux  passes  through  the  reed,  alternating  direction  with 
the  frequency  of  vibration.  Gaps  3 and  k are  taken  to  remain  constant  and 
equal.  The  alternating  flux  through  the  reed  Induces  an  ac  voltage  in  the 
coil  around  it  and  causes  a current  to  flow  in  a circuit  with  a load  resis- 
tance  R.  . To  optimize  the  design  parameters  we  need  an  analysis  of  the  r»ag- 
netic  circuit  of  the  transducer.  This  will  enable  us  to  calculate  the  energy 
dissipated  in  Ri  as  a function  of  the  gaps  between  poie  pieces  (fg),  the 
reed  thickness  (t),  the  reed  displacement  amplitude  (a),  and  the  frequency  of 
vibration  of  the  reed  (f). 

THE  HAGNETIC  CIRCUIT  EQUATIONS 

The  magnetic  circuit  can  be  replaced  by  the  analogous  electric  circuit 
as  shown  in  Fig.  2,  where  the  various  electrical  admittances  have  the  values 
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Fig.  2,  Schematic  Representation  of  the  Magnetic  Circuit  and  its  Electrical 
Analogue 
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of  the  relevant  Magnetic  permeances  (P)  of  the  reed  generator,  and  the  mag* 
ni tudes  of  the  currents  deduced  will  then  be  the  same  as  those  of  the  fluxes 
(4)  carried  in  the  corresponding  bravtches  of  the  magnetic  circuit.  The  clr- 
cult  analysis  is  straightforward: 


A a ii  ■ Ti  in 

VA  P,  ♦ P, 


(1) 


♦.  - 


‘ P.  +P3  A 


(2) 


4> 


2 P2  ♦ P3  A 


A . A 
It  I 


Substituting  (I),  (2),  and  (3)  in  (k)  yields 
A P»  P2  \ «t 

% “ vpj  ♦ p3  * p2  ♦ y \ 


(3) 


(M 


(5) 


PA  P.+P3  P2*P3 

CP,  » %)<>,  ^ %) 

A P,  ♦ P2  ♦ 2Pj 

Substitution  of  (7)  in  (S)  finally  gives  us: 
j»  * Pl  P2  \ I 

%"*t  <P.  ♦ P,  * P„  ~J  Pl4(P,  ♦ P,  ♦ 2P,) 


I 3 


r , , ;i*ur;r  -y 


(6) 


(7) 


(8) 


<p,  ♦ p,Wp2  ♦ Pj) 

, P,  Pj  \ , (P,  ♦ P,)(P2  ♦ P,)  a . . 

% m *%  W,  ♦ P3  • P2  ♦ f y ^P,  ♦ P3J(P2  ♦ Pj)  ♦ pJp,  ♦ P2  ♦ 2Pj)/  (9) 
P,(P2  ♦ PJ  - P2(P,  ♦ PJ 

**  " *t  <P|  ♦ Pj)(P2  ♦ ♦ P4(P,  ♦ P2  ♦ 2Pj)  (,0) 

P,(P,  - PJ 

**  " *t  (P,  ♦ P.)(P2  ♦ P3!  ♦ P^P,  ♦ 2P.  ♦ P2)  (,,) 
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THE  PERMEANCE  CALCULATIONS 


At  this  point  it  is  necessary  to  obtain  expressions  for  the  various 
permeances  P in  terms  of  the  reed  displacement  amplitude  (a)  and  the  config- 
urational parameters  which  can  be  conveniently  varied,  namely  the  two  gap 

lengths  l , and  2 , and  the  reed  thickness  t.  S*nce  i . is  always  kept 

go  gA 

equal  to  i „ we  denote  both  gap  lengths  by  2 , 

9“  0 

These  permeances  will  be  calculated  by  essentially  the  same  techniques 
used  in  Ref.  4.  These  consist  principally  of  approximating  the  actual  flux 
paths  by  circular  arcs  and  then  applying  the  permeance  formulae  for  these 
paths  as  listed  in  the  standard  text  by  Roters.^  Suitable  modifications  of 
these  general  formulae  must  sometimes  be  made  to  deal  with  the  peculiarities 
of  a particular  configuration.  The  actual  arrangements  for  the  prototype 
generator  in  both  the  isolated  and  in-shell  conditions  are  shown  in  Figs,  3 
and  4 respectively.  The  numerical  values  of  the  pertinent  configurational 
dimensions  are  somarized  in  Table  I, 

The  permeances  associated  with  all  of  the  flux  paths  which  do  not  pass 
through  the  reed  are  collectively  denoted  by  P4.  There  arc  five  Such  paths 
which  are  schematically  illustrated  in  Fig.  5*  and  are  calculated  immedi- 
ately below  along  with  the  other  permeances  illustrated  in  Fig,  5»  All 
calculated  permeances  are  in  inches. 

(0  Ppit  (Fig.  5-A) 


The  flux  path  here  is  in  the  form  of  a parallel  plate  magnetic  capacitor 
and  the  permeance  is  simply  given  by 


, n ate  Area  „ P _ (.7) (.365) 
P4  Plate  Separation  S ,395 


.647“ 


02) 


P_.  will  be  slightly  modified  by  the  presence  of  the  reed,  but  since 
the  reed  Is  very  thin  compared  to  the  separation  of  the  keeper  plates  this 
modification  is  very  small  and  will  be  neglected. 

(2)  (Fig.  5-B) 


This  is  a standard  permeance  path  for  which  the  formula  is 

.318  Y v1  +T  ) 03) 


where  y is  the  common  perimeter  length  of  the  keeper  edges  linked  by  the 
flux  path.  In  the  fluidic  generator  y 2,2  inches,  w is  the  edge  thickness 
.07"  and  / is  the  diameter  of  the  smaller  of  the  two  semicircular  boundaries 
of  the  flux  path  (.395").  Insertion  of  these  values  Into  Eg.  (13)  yields 

4.  H.A.  leupold,  F.  Rothwarf,  C.C.  Campagnoulo,  J.  Fine,  and  H.  Lee, 
"Magnetic  Circuit  Design  Studies  for  an  Inductive  Sensor,"  ECON 

Technical  Report-4158,  October  1973. 

5.  H.  Rotors,  Electromagnetic  Devices. (John  Wiley  fr  Sons,  Inc,,  1941). 
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Fig.  4.  Sectional  View  of  Fluidic  Generator 
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Py4  - .318  (2.2)  In  ^1  + 2la|Zl)  - (.318)  (2.2)  (.303)  - .212" 

(3)  (F?g,  5-C) 

This  Is  not  a standard  permeance  but  an  expression  for  It  can  He  derived 
as  follows:  Consider  for  the  time  being  only  that  portion  of  the  permeat.ee 
path  formed  by  the  annular  ring  of  width  W as  shown  In  Fig.  6.  As  we  shall 


Fig.  6.  Calculation  of  the  Permeance  P, 


MU 


see  a simple  extension  of  the  resulting  formula  can  then  be  used  to  determine 
the  remaining  permeance  consisting  of  the  semicircular  region  of  diameter  £, 
The  permeance  of  a path  such  as  the  annular  ring  Is  of  the  same  standard  form 
as  Eq.  (13) » l.e. 


P - .318  y Is^uji) 


m 


Similarly,  the  portion  of  the  permeance  within  the  annular  ring  of  width  w Is 


P - .318  V In  ♦ X1 ) 


9 


05) 


and  the  permeance  of  element  dw  is 
• 3 1 8y  7 dw 

JB  _ * 


2w 
I + 


However,  since  the  flux  tines  are  emanating  from  the  sides  of  a magnet, 
the  magnetomotive  force  driving  an  element  such  as  dw  is  only  a fraction  of 
the  full  HHF  of  the  magnet.  Since  this  MMF  is  in  direct  proportion  to  the 
distance  of  the  path  terminals  from  the  magnet  center,  this  fraction  is 
equal  to 

!♦» 

f-7 07] 

f.w 

and  the  effective  permeance  of  dP  is  reduced  by  the  factor  f,  viz.. 


’eff  * fdP  " -3I8V  ~i 


. '2 

l?  ■*  I 


I « . i » \ • t 

Kz  * yJ  <'*~j 


and  the  total  permeance  of  the  annular  ring  W is 


^ ! *. 


I8v  W 


"**  0 


To  obtain  the  total  effective  permeance  P^  one  needs  only  to  set 
l • 0 In  the  expression  for  P#^f  and  we  are  left  with 

^ ■ .3»8y  (2 

y is  just  the  perimeter  of  the  magnet  minus  the  length  of  the  inner  wall. 
The  latter  length  is  subtracted  because  it  is  rendered  largely  ineffective 
by  competition  from  the  flux  paths  of  permeance  P«^»  The  effective  y is 
somewhat  awkward  to  estimate  because  of  irregularities  such  as  the  screw 
slots  cut  into  the  magnet  walls,  hut  2,3"  seems  a reasonable  value  for  the 
two  magnets  together.  Accordingly,  our  value  for  P^  is 

• C»3 18)  (2,3)  - ,731" 


(«0  PU  (F{9*  S"0) 

This  again  is  a standard  permeance  of  the  form  used  in  Py*.  There 
is  one  such  path  for  the  outside  surface  of  each  pole  piece  so  we  must 
multiply  Eq.  (13)  by  two.  The  flux  lines  emanating  from  the  inner  surfaces 


go  to  the  reed  and  hence  contribute  to  Pjt  ?2,  and  Pj  rather  than  to  P^. 

For  y we  take  the  average  width,  (.1*75")  of  the  effective  portion  of  the 
trapezoidal  outer  surface  of  each  pole  piece.  2,  (.30")  is  twice  the  length 
of  the  reed  projection  beyond  the  gap  as  is  clear  from  Fig.  5-0.  As  .235" 
is  the  height  of  the  trapezoidal  pole  piece  plus  the  gap  from  the  pole  face 
to  reed,  it  also  follows  from  Fig.  5-0  that  w ■ (.235  - .150)  ■ .085".  in- 
serting these  values  in  Eq.  (13)  yields 

\l,  - 2<.JI8)(.MS)  In^l 

PLi.  ■ .302  In  1.57  - . 136" 

(5)  Pwi|  (Fig.  5-E) 


The  formula  for  a path  of  this  kind  is 


(21) 


and  since  we  have  four  such  paths  (2  for  each  pair  of  pole  pieces) 


(22) 


2 is  just  the  thickness  of  our  pole  pieces  which  is  .070"  and  A is  ~ n 

radians,  r is  the  distance  from  the  vertex  of  angle  0 to  the  larger  base 
of  the  trapezoidal  pole  piece.  Thus,  r2  is  equal  to  the  length  of  the  pole 

piece  edge  plus  one-half  the  gap  width  divided  by  the  sine  of  y radians. 

So  since  the  pole  piece  edge  is  ,2"  long  * 

I 

r2  • .2  ♦ CSC  y • .2  ♦ .577  (23) 

r(  is  equal  to  y4  esc  j plus  the  distance  along  the  pole  piece  edqe  from 

which  the  flux  lines  go  to  the  edge  of  the  reed.  The  latter  distance  Is 
just  one  half  the  reed  thickness  (t),  So  finally  we  have 


and 


r|  * 2 ♦ 'W  \ 


, . maaim  ;*  * -”7  S 

“*•  Jn  i ♦ ~j  I 

2 g 

.A  ♦ 1.10  a 

P . - • 1 3**  In 

t ♦ 1.16 


m 

(25) 

(26) 


It 


'-..V.'r-^TV  'K-,  p^,^nXVTCi*iA  /rr  W-5  * 


(6)  % 

We  now  have  expressions  for  all  of  the  components  of  the  "leakage 
permeance"  P^,  so  P^  is  given  by: 

P4  « PP4  + pyi*  * Pm4  + Pl4  + P«4 

A * 1.16  i 

Pu  - .647  + .212  ♦ .731  + .136  - .13*  In  t ♦ fjTi 

9 

A ♦ 1.16  ,* 

p*- 1-72  ♦ .»*  '"tttt rr  (27) 

9 


If  the  gap  is  varied  In  the  usual  manner,  i.e*  by  placing  pcrraal  loy 
shims  between  the  magnets  and  the  permalloy  keeper,  P^,  Pp^,  and  P^will 

be  slightly  affected,  but  since  the  shim  thicknesses  are  only  of  the  order 
of  10  or  20  thousandths  of  an  inch  the  change  in  these  quantities  will  be 
negligible. 


(7)  PG3  (p<9.  S-r) 


This  is  the  permeance  of  the  gap  between  the  reed  and  either  pole 
piece  at  the  pinned  end  of  the  reed.  The  length  of  this  gap  L is  given  by 


L • 


(28) 


Since  the  widths  of  the  pole  piece  (.313")  *nd  the  reed  .270"  are  not  the 
same,  we  take  an  average  (.291")  in  our  determination  of  the  gap's  cross- 
sectional  area.  The  pole  face  thickness  is  .07(7'  so 


Arc  2 (.070)  (.291) 
03  l 1,  • « 

(29) 

, . -<*08 
03  i,  - t 

(30) 

If  the  gap  sizes  are  varied  by  nwchining  the  pole  faces  the  constant 
in  this  formula  will  change  due  to  the  alteration  in  cross-sectional  area 
of  the  pole  faces* 

(8)  (fig.  S-C) 

P|*  is  the  permeance  of  the  path  extending  from  the  edges  of  the  reed 
to  the  pole  face  edges  parallel  to  the  former.  The  standard  expression  for 
this  type  of  permeance  path  is 

P • ,26y  (30 
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where  y is  the  edge  length  which  here  is  the  thickness  of  the  pole  piece 
(.070")  so 

P » (.26) (.070)  " .0182”  (32) 

Since  we  have  two  such  edges  the  total  permeance  is 

PE3  « 2P  * (2) (.0182)  « .0364"  (33) 

(9)  PF3  (rig.  5-H) 

These  are  the  paths  from  the  edges  of  the  pole  pieces  to  the  flat  por- 
tion of  the  reed.  These  are  standard  forms  whose  permeances  are  given  by 

P » .52y  - ,52(.29l)  (34) 

Since  there  are  two  paths  we  have 

PF3  « 2P  - 2(.52)(.29l)  - .303"  (35) 

(10)  PL3  (Fig.  5-1) 

These  permeances  involve  the  paths  between  the  walls  of  the  pole  pieces 
and  the  wide  faces  of  the  reed.  The  general  expressions  are  the  same  as 
those  for  P^  except  for  a factor  of  2 to  take  into  account  that  the  flux 
lines  here  are  quarter  circles  rather  than  semicircles  as  In  P^l*.  The  gap 
length  4 used  in  the  formula  must  be  taken  as  twice  the  separation  between 
the  reed  and  the  pole  piece,  i.e., 

*q  * * \ 

X - 2 <-5 — ) ■ *9  • * • 


The  quantity  w;  for  the  inner  flux  path  is  just  the  altitude  of  the 
trapezoidal  pole  piece  which  is  .15".  For  the  effective  y one  chooses  an 
average  of  the  reed  width  (.270")  and  the  width  of  the  trapezoid  midway 
*'*><een  its  bases  (.406").  This  results  in  a value  y ■ .338". 

The  outer  path  must  be  treated  slightly  differently  because  the  reed 
projects  'inly  .IS"  beyond  the  gap  edge.  Accordingly.  wQ  is  just  taken  as 

<'g  * ») 

wQ  - .15  - (36) 

and  the  averaqe  value  y is  also  slightly  altered. 

As  Is  seen  in  Fig.  7 the  ftux  lines  from  the  reed  extend  only  up  to 
line  AS  which  is  .15"  above  the  top  surface  of  the  reed.  The  average  width 
of  the  flux-emanating  region  of  the  pole  face  Is  then 


v - tLLjJU 

F 2 


(37) 


f«g.  7.  Determination  of  the  Width  yp  in  the  Formula  for  Permeance  P 


and 


V0‘ 


yF  * VR 


Where  yR  » .270  is  the  reed  width. 


i • t 

*B  - .)l)  * 2 {.1}  - v-2-j— ))  IM  30° 

AB  • .3*3  ■»  (.30  - (4g  - 0)  .577 

AB  • M 6 - .577  (I  - t) 

9 

M6  - .577  (/Q  - t)  * .313 

VF  T8 ■ .399  • .289  (>'  - t) 

.399  - .289  (lq  - t)  ♦ .270 

y„  • 1 ■ — 

y0  2 

y0  » .335  - . iaa  (;  - t) 

Substituting  these  values  It  the  general  expression  yields 
.l3  • 2<.3t*)(.3S8>  In  ^1  . 7^) 


♦ 2(.JI8):.335  - .!<*  (i  * t) ] In  i I ♦ 


♦ ^ .213  - .0916  (/  - o)  in  ( T^'t  ) (**) 

9 9 

00  Pw3  (Fig.  5-J) 

1 1 - . .arrow  side  face  of  the  reed  were  coplaner  with  the  narrow  face 

of  :hc  >ie  piece  as  indicated  by  the  dashed  line  in  Fig,  5-J,  the  expression 
for  permeance  Pwj  would  be  the  often  used  standard  form 


PwJ.  .318,  l,(l  *if) 


Although  it  is  possible  to  correct  for  the  fact  that  the  actual  flux 
lines  are  shorter  than  in  this  ideal  case,  the  correction  would  be  very  small, 
and  since  the  total  permeance  is  itself  so  small  as  to  barely  affect  the 
calculations  with  one  significant  figure,  this  correction  is  completely 


negligible  and  will  not  be  made  here.  So  w is  taken  as  j,  4 • j 


lq  - t 


cos  30®  • 


y * .070,  and  we  precede  the  expression  with  the  factor  2 because  there  are 
two  paths  of  this  type 


Pw3  - 2 ( .3 18) (.070)  In  (l  ♦ -y”  ^ ) 
* -0M.5  In  V1  * rTT  ) 


(12)  8C,  W8  PM 


These  permeances  are  of  the  same  form  as  PC3  with  allowance  made  for  the 
spacing  changes  due  to  the  displacement  of  the  reed.  Mien  the  reed  Is  at 
equilibrium,  that  is,  at  zero  displacement,  PqI  * *G2  ■ *C3«  Since  we  are 
interested  in  the  max I nun  and  mfnlmum  values  of  Pqi  and  Pq2,  we  write  them 
In  terms  of  the  displacement  amplitude  (a)  of  the  reed. 


1 - t * 2a 

9 


pnax  #ma 

ul  G2 


pm In  pain  ,0408 

Gl  rC2  4 • t ♦ 2a 
9 


V*,l>  <'*-  •■'i-K-Vfr^:  « 


where  P 


Gl 


is  a maximum  whe:  P^ 


is  a minimun  and  vice  versa. 


03) 


PE!  ‘ PE3  ’ •°36'*' 

<"*>  fJl 

PFI  * PFJ  ’ •*»" 

05)  pw)  — Pw! 


(52) 


(53) 


PW|  is  of  the  same  form  as  Pwj  when  the  reed  is  in  equilibrium.  When 
Pw|  is  a maximum,  1 ■ (ig  - t)  must  be  replaced  by  (lq  - t - 2a)  and  by 
(*g  * t ♦ 2a)  when  Pwj  is  a minimum.  The  same  holds  For  Pw2.  Thus  we  have 

(5M 

(55) 

(16)  PL|  and  PL2 


p*ax 

w 1 


max, 

w2 

min 


.0445  In  [ 


I ♦ 


1.73t 


„min 

P , ■ P . 

w I w2 


L - t - 2a  I 
9 

,,  ».73»  n 

•0,A5  l"(l>r-7i,i 

9 


The  term  in  Pjj  arising  from  the  Inner  flux  path  is  the  same  as  the  cor- 
responding term  in  PLj  except  that  again  tha  quantity  (/g  - t)  is  replaced 

by  0*g  - t - 2a)  or  (/g  - t + 2a)  to  obtain  pJJ*  and  P^j",  respectively. 


The  term  arising  from  the  outer  flux  path  is  similarly  derived  from  the 
corresponding  term  in  P|_j  by  the  same  substitution  but  an  additional  altera- 
tion is  necessary.  The  portion  of  the  reed  projecting  beyond  the  limits  of 
the  gap  is  triangular  rather  than  rectangular  as  is  the  projection  at  the 
anchored  end  of  the  reed.  For  this  reason  the  flux  lines  emanating  from  the 
walls  of  the  pole  piece  have  only  one  half  the  reed  area  in  which  to  termin- 
ate. To  take  this  into  account  it  seems  reasonable  to  substitute  in  the 
expression  for  yg,  one  half  of  the  reed  width  for  the  full  reed  width.  Thus, 
in  Eq.  (42)  the  value  of  .270  is  replaced  by  .135.  After  all  these  adjust' 
ments  are  made  the  expression  for  P^  and  P^2  becomes 


max 


max 


pmin  plain 
’ll  " *12 


.215  In  ♦ 


-J1 


f - t * 2a  ) 
9 


♦ '.170  - .0916  (/  - t * 2a)  1 In  7 «j°x  , 

g Zg  - t » 2a 


(56) 


There  Are  other  permeances  associated  with  flux  paths  going  from  comer 
to  corner,  edqe  to  comer  and  face  to  comer.  In  this  problem,  these  are  not 
of  standard  form  and  are  very  difficult  to  estimate.  They  tend  however,  to 
be  small  compared  to  most  of  the  calculated  permeances,  and  the  labor  of  any 
attempt  to  estimate  them  would  not  be  justified  by  the  precision  attainable 
by  our  computational  techniques* 
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CALCULATIONS  OF  GAP  FLUX  DENSITIES  IN  THE  ABSENCE  OF  A REED 

We  now  have  expressions  for  ell  of  the  permeances  relevant  to  the  de- 
sired power  calculations.  Before  a sample  calculation  for  particular  values 
of  4g»  t and  a are  made,  it  would  be  instructive  to  calculate  the  field  In 
the  gaps  when  the  reed  is  not  present.  Since  a measured  field  value  exists 
for  a gap  width  of  .060,  a comparison  with  the  result  of  our  calculation 
would  give  an  indication  as  to  the  general  efficacy  of  the  computational 
methods  used.  To  do  this  we  need  the  permeance  Pg4  which  is  given  by: 

, . (2)  (Are.  Pol.  F,c«s)  _ (2H.02I9)  . (57) 

rG4  Gap  Length  Too 

We  now  also  have  the  permeances  Pfi*  between  the  corresponding  edges  of  the 
pole  faces  and  Pc4  between  the  corresponding  corners  on  either  side  of  the 
gap.  P|^  is  given  by  the  usual  formula  P £4  ■ 2(.26)y,  where  y is  the  perim- 
eter of  a pole  face  and  the  factor  of  2 enters  because  there  are  two  gaps. 

Pc4  is  given  by  (8)(.077)£g,  where  the  factor  8 occurs  because  there  are 
four  corner-to-corner  permeances  for  each  of  the  two  gaps.  Therefore,  Pc4 
is  given  by 

P^  - (8) (.077) (.060)  • .0369-  (58) 

PE4  - (.52)  (.766)  - .398.  Pp4,  Py4  and  Pj$4  are  the  same  as  when  the  reed 
is  present.  We,  therefore,  need  only  to  find  the  new  values  of  fyi*  and  P^4* 
Pgj*  is  given  by  Eg.  (26)  with  t set  equal  to  zero.  Thus 

•U  ■ •'* '» **■ ■ •'» ,B  6-7i 

P^  - (.13^(1.91)  • .256"  (59) 

The  contributions  of  the  two  outer  surfaces  to  P^i,  are  given  by  the  same 

formula  as  when  the  reed  is  present,  but  now  1 - ig  « .06  and 

y ■ j " *606.  these  contributions 

P°l*  ■ 2(. 318) (.406)  In  ^1  ♦ 2ig|I  j . .258  In  7.67 

P°lk  - (.258)  (2.04)  . .526"  (60) 

In  the  absence  of  the  reed  the  inner  surfaces  also  contribute.  The 
longest  semicircular  arc  of  this  path,  however,  must  be  no  longer  than  the 
distance  d as  shown  in  Fig.  8,  so 

*'•)'*  (61) 

d £ 

" " n " 2 (62) 
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Ftg.  8.  Determination  of  the  Inner  Flux  Paths  Comprising  Permeance  P 


f 


w • .126  - .030  • .0958’ 
and 

5 . Jl4Ul.-t.Mj.mP>  »"  30°  . .J,J  . .o:5  - .366" 

Thus 

• 2(.3I8)(.J68)  In  vl  ♦ ) - .*$<♦  In  (4.03)  » .325 

\k  m PLh  * PU  " *325  4 *526  " •:5,‘ 

Our  total  penmeanee  P^  in  the  absence  of  a reed  is  then 

\ ■ ftk  * '*  * * \t  ♦ 'tk  * fdk 

p,.  • .730  ♦ .647  ♦ .712  ♦ .731  ♦ .256  ♦ .851  ♦ .398  ♦ .037 


mag* 


In  this  case  P.  Is  our  total  penaeance  P and  from  It  we  can  obtain  our 
load  line  slope 

<m> 

Where  L - .395  In  and  A - .337  In2.  And  from  the  intersection  of  the 
Alnico  5 demagnetization  curve  with  the  load  line  B/H  ■^+.52  in  Fig.  9 we 
obtain 

B ■ 2670  gauss  (68) 

mag  3 

Since  the  magnetic  configuration  Is  of  irregular  shape,  B is  not  perfectly 
uniform  over  the  entire  magnet,  therefore,  B^g  represents  an  effective 
average  flux  density.  So  the  total  flux  output  of  the  magnet  is 


- BA  - (2670) (2.1 7)  - 5790  Maxwells 
t AM 9 M9 

The  flux  in  one  gap  4>g  is  given  by 

p 

♦c  • 2^  *t  ■ 5790  • w 

and 

Bc  " ^ ^ " 3880  9*u*$ 


(69) 

(70) 
(70 


This  value  is  about  29%  higher  than  the  measured  value  of  3000  gauss 
which  is  satisfactory  agreement  considering  the  various  approximations  made 
In  the  course  of  the  calculations.  One  must  also  keep  in  mind  that  such 
comparisons  should  not  be  taken  too  seriously.  The  calculated  Bg  represents 
an  average  flux  density  obtained  by  dividing  the  total  flux  emanating  from  a 
pole  face  by  the  area  of  that  pole  face.  The  measured  Bg  can  represent 
either  a much  more  localized  value  of  Bg  or  on  average  that  Includes  fluxes 
not  even  associated  with  the  gap,  depending  on  the  type,  positioning  and  size 
of  the  probe  used.  The  wider  the  gap  compared  to  the  pole  face  dimensions, 
the  more  difficult  It  is  to  obtain  a meaningful  measured  value  of  Bg.  Con- 
sidering these  limitations,  we  conclude  that  the  general  computational  proce- 
dure yields  reasonable  agreement  with  the  measured  results,  and  that  Its 
judicious  application  to  the  determination  of  expected  power  outputs  should 
not  lead  us  too  far  astray, 

SAMPLE  CALCULATION  FOR  POWER  OUTPUT  IN  THE  LOAD  RESISTANCE  I ^ 

We  will  now  perform  a sample  calculation  for  the  most  efficient  config- 
uration measured  thus  far.  The  relevant  parameters  for  this  arrangement  are: 

- .070'* , t - .032",  a - .01?*. 

9 


We  will  first  calculate  P4.  We  need  only  the  components,  P„j,  and  P.l. 
Using  Eq.  (26)  we  obtain  1 
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V " J34  ,n  4‘26  " •|3^(K^5)  - .194"  (72) 

The  appropriate  Pl4  has  already  been  calculated  on  page  II  and  found  to 
be  ,*36.  The  other  components  of  P4  are  tr.a  sane  as  were  used  to  calculate 
the  gap  fields  except  tha^.  now  there  is  no  Pc4,  Pfi*  or  Pc^.  Thus 

P4  " Pp4  + Py4  * PH4  + PL4  * Pw4  ™ 

P4  - .647  + .212  + .731  ♦ .136  * .'94 

“ 1.92"  (74) 

p , has  five  components 


PG3  • T070°^032  - <W 

PE3  • .0364* * (76) 

PFJ  - -i03"  (77) 

Pu--J,s  ♦7J?t} 

3 

* (.213  - .0916  (l9  - 0)  In  ) (78) 

- .215  (•  ♦ ) ♦ (.213  - .09'6(.038»  In  (^  ) 

pu  - .215  In  8.90  * (.213  - .0035)  In  7.90  • (.2151(2.19)  + (2.09X2.07) 
PL,  • .1*70*  ♦ ,1*32" 

\j  ■ -S02"  (79) 

P.)  - .01*5  In  0 ♦ } * -0«.5  I.  0 ♦ ) (8J)) 

Pw|  * ,n  " (*0445) (.900)  • .0400*  (81) 

Pj  • 1.07  ♦ .0564  ♦ .303  ♦ .902  ♦ .0400 

Pj'-T.iv  (82) 


>n<£ 


Pj  and  P2,  when  the  reed  is  at  its  equilibrium  position,  are  obtained  by 
substituting  a * 0 in  the  expression  for  the  components  of  P|  and  P2  and 
then  adding  these  components.  This  procedure  yields  P|  ■ P2  “ 2.27 We 


need  to  find  p|,  that  is,  the  total  permeance  at  equilibrium  to  locate  the 
point  on  the  Alnico  5 demagnetization  curve  at  which  the  recoil  line  begins. 
(Point  E in  Fig.  9) 


E Pl 

Pt  ?k  * 2 + T " 4,23 


| pt  “ ‘ i»J7)(H.23)  - - M5 


6^  =2950  gauss  (from  Fig.  9) 


H 


2950 


05 


* - 596  Oe 


(83) 

(84) 

(85) 

(86) 


The  slope  of  recoil  lines  for  Alnico  5 is  4.3  so  the  equation  of  the 
recoil  line  is 


B - 4.3H  ♦ B, 


(87) 


2950  * (4.3)(-  596)  ♦ 60 
2950  ♦ (4.3) (596)  - Bq 
2950  ♦ 2560-  #q 
>510  • B„ 


(88) 


So  we  have  for  our  recoil  line 
B » 4.3M  ♦ 55*0  gauss 


(89) 


We  now  need  on!/  Pj  and  Pj  to  complete  our  analysis.  Using  Eqs.  (50) 
and  (51)  we  have 


.max  _ „max 


61 


a — . .0408  .0408  . .0408 

62  £ - t - 2a  .070  - .032  - 2(.OI7)  toB* 

9 


Dmax  pmax 

«l  G2 


<0.2" 

.0408 


a&in  m j/*in  — .0408. — _,Q408  — » caa>< 

P62  FGI  1 - t ♦ 2a  .070  - .032  ♦ .034  .072  *5 


(90) 

(91) 

(92) 


From  (33)  and  (35)  respectively,  we  have 


E I 


- P, 


£2 


.0364* 


(93) 
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mm 


**  • - v-*  wa*  5 


p « P ■ .103 " 
rF\  F2  33 


(94) 


Fro*  (54) 


*1 


,«ln 


w2 


rj .[tc  tn  ( 1 -f  N)  m obUS  in  ( 1 + » — ^ lZ<LlXxfi22JL  ■ ) 

.0445  In  [l  ♦ t — * 2a  ; - .W  in  ♦ .070  - .032 “ (.01 7)  ) 


P™x  - .Okks  In  (l  ♦ ^55^  ) - (.0445)  In  14.85 


*£?*  " (’^S)  (2.70)-  .120 1 


p*2n  " 'okks  ,n  ♦ 0,i7o^  ) “ MS  ,n  1,77 


.min 


P"-"  . (.0445) (,57*)  - .0254 


From  (56) 

P~x  - .215  In  [\  ♦ «•  [.170 


(.0916) (.004)1  In 


P[7*  - .215  In  76  ♦ .170  In  75  - (.2I5)(4.33)  ♦ (.I70)(4.32) 


Pft*  * .931  ♦ .734  - 1.66 


LI 


Similarly 


P*J"  • .215  ^ ♦ ^072  ) * M70  - (.09»6)(.072)]  In  ) 


Pj*n  - .215  In  5.17  ♦ .163  In  4.17  ■ .215(1.64)  ❖ (.163)0.43) 


Pjj"  • .353  ♦ .233  ■ .586" 


.•In 


Finally,  we  have  for  P(  and  P2 

P^**  • 10.2  ♦ .036  ♦ .j03  ♦ .no  * 1.663 
P?**  • 12.32* 
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P"""  * .566  + .036  + .303  + .025  + .586 
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p"«*  c , 92  «.  U±*21  tJjJ&llhSl  1 2f3?) 
t 0 12.32  + 2.35  + 1.52  + 2.35 


p***  . , 92  + UMlUhZll  m k gg,, 

Kt  18.5  *** 


B N, 

Vh  J 


max 


J P^X  - - ( 1 * 1 7)  (4.98)  - - 5.83 


From  (89) 


Bmax  * 4.3H  + 5510 


„max  5510  5510  _ 

B * *gmu.*  * im  3*70  gauss 


1 + 


1.73 


5.83 


4P8*  - B1'”"'  A 


^ ‘magnet  * foWHaje)  » 6880  Maxwells 


Substituting  all  of  our  calculated  permeance  values  and 
in  Eq.  (II)  yields  for  the  amount  of  flux  carried  by  the  reed  1 


«T  * 6880  ( W*. 7)  13^7)  *+2(  1.92? 08.5) 

*™ax  . (^SQ5)C|ff3£)  . 1890  Maxwells 

The  Induced  voltage  amplitude  is  then  given  by 


d/n#x 

."tax  _ N k it 
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8 


where  is  the  angular  frequency  of  vibration  of  the  reedt  which  in 
caso  is  2^(1500)  and  N ■ (1500)  is  the  number  of  turns  in  the  coil. 


vmax  . 1LS00 . 267  volts 
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Attached  to  the  coil  is  a 600-4}  load  resistor  and  the  resistance  of 
the  coil  is  600  ft  as  well.  The  inductance  of  the  coil  with  a .01?'  reed 
is  about  23  mhenrys.  To  find  the  power  dissipated  in  the  load,  we  need  the 
total  impedance  Z of  the  circuit. 

Z2  « A2  ♦ X?  - {1 200) 2 ♦ X?  (109) 

L L 


The  inductive  reactance  X.  is  just 


«L  ■ 2T.fi,  ■ . 2I7  ti 

(110) 

Z2  « ( 1200) 2 ♦ (21 7)2  - 1.49  x I06 

(Hi) 

Z * 1220  ft 

(112) 
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The  root  mean  square  current  T is  given  by  I « 


max 


max 


Jl 


1/2 


and  the 


power  dissipated  in  the  load  is  then 


l2A. 


«* 


600 


2) (1220) (1220 


!)  ■ 


k watts 
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The  actual  power  measured,  however,  was  only  5 watts.  The  foregoing 
calculations  were  made  on  the  assumption  that  the  reed  had  infinite  permeance. 
The  electrical  steel  of  which  the  reed  is  composed  saturates  at  19,900  gauss 
and  its  cross-section  area  is 

A • (.2?  x .032  x 2.51*  x 2,54)  cm2  (114) 


Therefore,  the  maximum  amount  of  flux  it  can  carry  Is 

- ( 19*900) (.2 7) (.032) (2 .5**) (2. 5M  • 1110  Maxwells  (115) 

and  if,  as  was  shown,  1890  Maxwells  in  the  reed  wot: Id  yield  14,4  watts  then 
the  maximum  allowable  flux  of  1110  Maxwells  would  give  us 

v 2 

Ps.v  jj~”“  j llt.l*  « lt,97  watts  which  is  almost  exactly  the  power  observed. 


It,  therefore,  appears  that  to  utilize  the  full  flux  output  of  the  mag- 
nets used,  the  reed  thickness  must  be  Increased  by  about  70&.  Whether  this  can 
be  done  without  drastically  affecting  the  vibrational  response  to  the  oscilla- 
tor is  a mechanical  problem  beyond  the  scope  of  this  study.*  It  is  clear, 
however,  that  increased  output  with  the  present  configuration  depends  largely 
upon  increasing  the  flux  carrying  capability  of  the  reed.  If,  on  the  other 


3*  Carl  J.  Campagnoulo  and  Richard  N.  Gottron,  "The  Fluidic  Generator:  A 
New  electrical  Power  Source,"  24th  Annual  Proceedings  Power  Sources 
Symposium,  May  19*21,  1970. 


hand,  the  present  output  of  5 watts  ts  satisfactory,  it  can  be  attained  with 
considerably  less  magnetic  material  since  with  the  present  arrangement  7 0% 
•ore  flux  Is  available  than  can  be  profitably  used.  Of  course,  if  the  reed 
thickness  is  increased,  the  gaps  must  be  lengthened  by  an  appropriate  amount 
to  keep  the  load  line  of  the  magnets  the  same  and  to  maintain  reasonable 
space  in  which  the  reed  can  oscillate.  Also  the  expected  power  output  would 
be  altered  slightly  since  the  Inductance  of  the  coll  is  a function  of  reed 
thickness.  This  effect,  however,  would  not  be  significant. 


Parts  of  the  keeper  may  also  be  saturated  under  the  given  operating 
conditions.  The  cross  sections  which  would  saturate  first  are  the  ones  in- 
dicated in  Fig.  3*  Under  the  given  conditions,  the  flux  that  must  be 
carried  by  the  narrowest  sections  of  the  keeper  Is  given  by 


and 


c*  • v * v 

2 

pmax 

t 

6880 

kj8  - .731  - .212 

2 

A.98 

6880 

2 

• 2790  Maxwells 

11  • I • • ly>M0 


(116) 

(M7) 

(118) 


Since  permalloy  saturates  at  16,000  gauss,  the  keeper  cross  section 
need  be  increased  by  only  7%  to  take  full  advantage  of  the  available  flux 
at  the  given  reed  displacement.  This  change  alone,  of  course,  would  do  no 
good  unless  the  reed  thickness  were  doubled  as  well,  as  we  have  already 
shown. 


EXPERIMENTAL  CHECK  OF  CALCULATIONS 


To  test  the  validity  of  the  foregoing  calculations,  a series  of  meas- 
urements was  made  to  determine  the  maximum  voltage  output  as  a function  of 
reed  displacement  under  conditions  for  which  there  is  no  magnetic  satura- 
tion. The  reed  displacement  was  measured  by  "stopping"  the  vibrating  reed 
at  its  maxi max  displacement  (a)  with  a strobe  light  and  then  measuring  (a) 
with  a telemicroscope  equipped  with  a calibrated  micrometer  screw. 

The  results  of  this  experiment  are  shewn  in  Fig.  10  where  they  are 
compared  with  the  theoretical  curve  plotted  from  results  of  the  computer 
calculations  listed  In  Appendix  II,  As  Is  apparent  from  Fig,  10,  agreement 
Is  excellent  considering  the  many  approximations  used  In  the  calculations 
and  that  (a)  can  be  measured  with  a precision  of  no  better  than  ten  to 
fifteen  percent. 
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fiq.  10.  Comparison  of  Experimental  Points  (dots) 
with  Theoretical  EMF-Displacement  Curve 
(sol  id  curve) 


LOSSES  IN  THE  RICO 

It  is  of  interest  to  know  how  much  electrical  energy  It  dissipated  In 
the  reed  through  eddy  current  losses  and  hysteresis.  The  expression  for  the 
eddy  current  power  dissipation  In  watts  (see  Appendix  I)  is 


-Vl’-5;. 


6c  10 


TT 
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6o  10 


16 


019) 


where  n is  number  of  laminations,  f is  the  frequency,  II,  is  the  flux  density 
amplitude  in  gauss,  * is  the  lamination  thickness,  < and  w are  the  reed 
length  and  width  respectively,  V is  the  lamination  volume,  and  ; Is  the 
resistivity  In  ohm-cm.  The  reed  dimensions  are  all  measured  In  centimeters 
and  the  power  is  In  watts.  Since  the  reed  is  composed  of  two  laminations 
each  of  the  thickness  * • .0k06  cm,  we  multiply  (119)  by  two  and  insert 
.0M>6  for  ' at  well  as  the  appropriate  values  of  the  other  parameters  under 
saturation  conditions. 
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The  expression  for  hysteresis  power  loss  in  the  reed  is  given  by 


B«  J.6 


B .1.6 


? m % * 7 * 

p«'Tf  ^ : £"‘  ' T,f  v 


(121) 


where  V • .0706  cc.  ’ is  a constant  depending  on  the  material  and  is  .002 
for  our  reed.  Halting  the  proper  insertions  gives  us 

v 1.6 


PH  « (.002)  (1500)0600)  v &&  * (.0706) 


10' 


o • ( .002) 0 500) ( 1 500) ( 1 200) ( .0706) 
I010 


* .COnoi.  which  it  netiliqib'e. 
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The  eddy  current  losses,  however,  could  become  quite  significant  if  the 
reed  were  not  properly  laminated.  For  example,  an  unlaminated  reed  of  the 
same  thickness  as  the  laminated  one  would  have  four  times  the  eddy  current 
loss.  It  is.  therefore,  important  that  if  the  reed  thickness  is  increased 
as  recommended,  it  be  done  by  adding  laminations,  A solid  un laminated  reed 
of  twice  the  thickness  of  the  one  presently  used  would  have  a power  loss  of 
about  i watts,  whereas  a fourfold  lamination  of  such  a reed  would  reduce 
the  loss  to  about  .68  watts. 

LOSSES  IN  THE  KEEPER 


Eddy  current  losses  in  the  keeper  come  about  because  of  fluctuations 
in  the  flux  density  as  the  magnet  load  line  changes  with  the  vibration  of 
the  reed.  If  the  keeper  thickness  is  increased  by  the  required  7%»  or 
better  still,  widened  at  the  narrow  points  which  saturate  by  7%,  the  maxi* 
mum  Mux  density  would  be  no  more  than  16,000  qauss.  The  operating  condi- 
tions would  then  be  optimized  and  It  would  be  useful  to  have  some  estimate 
as  to  the  edoy  current  losses.  Such  an  estimate  would  be  very  rough  since 
it  is  very  difficult  to  determine  the  flux  pattern  in  the  keeper.  We 
already  know  that  when  the  reed  is  at  its  greatest  displacement  the  flux 
carried  by  the  narrowest  portions  of  the  keeper  would  be  16,000  gauss.  We 
also  need  the  flux  density  when  the  reed  is  in  Its  equilibrium  position. 
This  is  obtained  by  substituting  the  appropriate  equilibrium  values  of  the 
parameters  in  Eqs,  (116  • 118),  In  this  fashion  we  obtain 
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The  cross-section  dimensions  of  the  keeper  are  .91**  cm  and  .178  cm, 
so 


>K  ■ t.9l4H”7»)  * X5’300  ,aU“ 

Therefore,  the  field  at  the  narrowest  portion  of  the  keeper  fluctuates  from 

15*300  gauss  to  16,000  gauss  with  an  amplitude  of  — i — * 350  gauss.  The 
fluctuation  in  other  parts  of  the  keeper  would  be  somewhat  smaller,  but 
330  gauss  can  be  used  to  get  an  upper  limit  for  the  expected  losses.  To 
this  end  we  use  Eq.  (119)  with  the  appropriate  constants  for  the  keeper 
whose  volume  is  2.****  cc  and  whose  thickness  would  be  .178  cm.  Further,  the 

fluctuation  frequency  of  BE  is  double  that  of  the  reed  since  BE  is  a maximum 

at  both  ends  of  the  reed's  vibrational  cycle.  Using  these  values  in  Eq. 
(119)  yields: 

p , (3.1**)  (3.  !*♦)  (3000)  (3000)  (350  )(350)(  .178 )( .178)  (2  .*♦*♦)  IQ6 
E (6)  (**5)  to’6 

p^  ■ .31  watts  (12**) 


which  is  about  2\  of  the  expected  optimum  output.  This  loss  of  electrical 
energy  can  be  avoided  by  a double  lamination  of  the  keeper  to  cut  the  loss 
to  U.0M  watts  or  even  a quadruple  lamination  which  would  result  in  the 
low  eddy  current  loss  of  0.0?  watts.  Since  these  numbers  represent 
only  upper  limits,  the  actual  values  would  probably  be  much  smaller.  The 
use  of  high-resistivity,  high-permeability  composite  castings  for  the 
keepers  would  probably  eliminate  the  eddy  current  losses  entirely. 


The  hysteresis  loss  in  the  keeper  cannot  be  obtained  directly  frem 
Eq.  (121)  because  the  field  in  the  keeper  varies  only  over  positive  values 
of  B,  while  (121)  applies  only  to  complete  cycles  about  the  principal 
hysteresis  loops.  We  can,  however,  get  some  idea  as  to  the  order  of  magni- 
tude of  the  loss  by  substituting  our  amplitude  350  gauss  Into  (121).  Using 
the  value  11  ■ .0001  for  permalloy  we  have 


PH 

PH 


(.0001)  (3000)  (3000)  ; tt*1*) 

k,X  x I0-6  watts 


(125) 
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which  as  in  the  case  of  the  reed  is  negligible. 


SUMMARY  AND  RECOMMENDATIONS 

(I)  Under  the  present  conditions  of  operation  the  fluidic  generator 
is  producing  only  about  3516  of  its  potential  power  output.  Ooubling  the 
cross-section  of  the  reed  by  adding  lamination  or  increasing  its  width 
would  augment  its  flux  carrying  capability  to  where  the  potential  power 
output  could  be  fully  realized  with  only  negligible  eddy  current  and 
hysteresis  losses  in  the  reed.  Appropriate  adjustments  in  gap  length  would 
also  be  necessary  to  prevent  saturation  and  to  allow  sufficient  latitude 
for  vibration. 
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(2)  Keeper  hysteresis  losses  are  also  negligible  but  eddy  current 
losses  can  be  appreciable  unless  the  keeper  is  laminated  or  made  of  a suit- 
able composite  material.  Further,  the  narrowest  part  of  the  keeper  cross 
section  must  be  increased  by  about  7%,  to  avoid  saturation  with  the  given 
reed  displacement  amplitude  of  .0 1 71 ' . 

(3)  To  facilitate  design  optimization  of  future  generators  a computer 
study  was  done  to  obtain  a series  of  design  matrices  which  list  voltage 
amplitude  as  a function  of  gap  length  (/g),  reed  thickness  (t),  and  reed 
displacement  amplitude  (a).  The  computer  program  uses  the  mathematical 
procedure  of  the  sample  calculation  made  in  this  report  and  it  is  written 
in  Fortran  IV  for  the  Burroughs  5500  Computer.  The  program  is  reproduced 
in  Appendix  II  together  with  a suninary  of  the  equations  used  in  the  cal- 
culations. 


The  computed  voltage  amplitudes  are  listed  in  Table  II  in  Appendix 
II.  Al I of  the  listed  voltages  are  calculated  for  a 1500  turn  coll,  a reed 
width  of  .2 T',  a reed  vibrational  frequency  of  1500  hertz, and  the  standard 
magnet  and  keeper  dimensions  used  in  the  sample  calculation  and  summarized 
in  Table  I.  The  assumption  is  also  made  that  the  keeper  does  not  saturate 
for  any  of  the  parameter  choices  listed.  This  assumpt ioi\  as  we  have  shown 
in  the  sample,  is  not  valid  for  the  larger  values  of  (a)  unless  the  keeper 
thickness  is  increased  by  an  appropriate  amount.  This  thickening  of  the 
keeper  would,  of  course,  affect  the  value  of  Py4  and  hence  the  permeance 
calculations.  The  effect  on  the  calculated  voltage  output,  however,  is  not 
likely  to  be  significant  because  in  the  sample  calculation  Pyt*  is  only  of 
the  order  of  10%  of  Pi*,  and  P4  In  turn  constitutes  only  about  kO%  of  P*. 
Furthermore,  Pyq  itself  is  relatively  insensitive  to  keeper  thickness  (W) 
since  W appears  only  in  a small  correction  to  t he  logarithmic  factor  in 
Eg.  (13)*  Hence,  thickening  of  the  keeper  to  required  values  will  probably 
not  affect  the  calculated  voltages  by  more  than  the  order  of  a percent  or 
two  even  in  the  worst  cases. 

The  calculations  also  assume  Infinite  reed  permeability  and  the 
values  not  actually  attainable  with  the  reed  material  used  are  marked  with 
an  asterisk.  Thus  the  design  can  easily  be  optimized  by  choosing  the  param- 
eters yielding  the  largest  non-asterisked  voltage  in  the  Tables.  This,  of 
course,  is  always  subject  to  the  condition  thet  the  reed  of  requisite  thick- 
ness can  be  driven  with  the  steted  amplitude  and  frequency.  Perhaps  tMs 
can  be  assured  in  many  cases  by  making  appropriate  adjustments  of  the  reed 
stiffness  throuqh  variation  of  the  reed  material  and  dimensions  at  the  point 
where  it  is  anchored  or  by  widening  the  reed  rether  than  thickening  it* 


The  Tables  can  be  used  to  obtain  the  power  output  in  eny  load  re- 
sistance A(_  from  Eq.  (126),  viz. 


V'  A. 
max  L 
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where  l Is  the  total  circuit  impedance  obtained  from  Eqs,  (109)  and  (110), 
The  coil  inductance  is  a function  of  reed  thickness  and  the  emplricat  curve 
of  the  relationship  is  shown  in  Fig.  II.  As  is  evident  from  (III),  however, 
the  reactive  portion  of  the  impedance  is  only  about  a percent  and  a half  of 
the  total, and  the  variation  of  reactance  with  reed  thickness  will  not  have 
a significant  effect  upon  power  output  at  the  frequency  and  resistance  used, 
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REED  THICKNESS  (T)  INCHES 


Fig.  n.  Inductance  or  Sensing  Coil  as  e Function  of  Reed  Thickness. 

(4)  The  present  configuration  of  the  Magnetic  circuit  results  in  rela- 
tively low-sloped  load  lines  which  intersect  the  Atnico  5 demagnetization 
curve  well  below  its  "kne#'  (Fig.  9),  Note  that  if  Co^Sn  were  used  in  place 
of  Alnlco  5 <*uch  higher  flux  densities  would  result  for  any  given  load  line 
with  a concomitant  increase  In  power  output.  This,  of  course.  Is  always 
assuming  that  neither  reed  nor  keeper  are  saturated.  It  probably  would  be 
most  practical  to  keep  the  power  output  at  about  10-12  watts  with  the  use  of 
much  less  magnetic  material  through  substitution  of  C t£$n,  A more  compact 
generator  operating  at  lower  toad  lines  with  the  sane  power  output  could 
probably  be  designed  using  the  latter  material. 

(5)  A composite  material  consisting  of  permalloy  particles  densely 
packed  in  a nonconducting  binder  night  be  desirable  for  the  keeper.  Pro- 
viding a sufficiently  high  packing  density  Is  obtainable,  eddy  current 
losses  could  be  virtually  eliminated  without  too  great  a loss  in  flux- 
carrying capability.  Such  a composite  could  also  be  press  molded  into  the 
desired  shape  thus  improving  the  efficiency  of  fabrication. 
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APPENDIX  I 

DERIVATION  OF  THE  FORMULA  FOR  ENERGY  LOSS  THROUGH  EDDY  CURRENTS 

Consider  a sheet  of  metal  such  as  the  one  pictured  In  Fig.  12.  The 
length  i is  much  larger  than  the  thickness  r and  the  width  w is  arbitrary. 

A sinusoidal  field  B is  acting  in  the  y direction  as  shown.  The  varying 
field  B will  induce  EMFs  which  will  give  rise  to  circulating  or  eddy  currents. 
We  wish  to  find  the  p<a*er  expended  by  these  currents.  We  proceed  by  consider* 
ing  a current  loop  Indicated  by  the  arrows  in  Fig.  12.  The  EMF  around  this 
loop  is  then  gi*en  by 


d»  _ AdS  _ 2/x  dB 
I08dt  lO^dt  I?  dt 
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if  we  neglect  the  small  voltage  drops  across  the  width  of  the  loop  at 
the  ends  we  can  write  for  the  voltage  drop  per  unit  length  in  the  z direction 
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The  current  density  i is  then  obtained  from  the  relation 


where  o is  the  resistivity  of  the  material.  Substituting  (129)  in  (130) 
yie»ds 


030) 


"To*' 


cos  u>t 


030 


The  instantaneous  power,  dP,  dissipated  in  the  unit  laminar  element  dx  is 
then 
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and  the  total  power  loss,  P,  of  a unit  length  of  sheet  is 
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The  energy  loss  per  cycle  is  then  given  by 
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,3  p 


— I cos2  (u>  t)  (<a  d t ) 
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The  average  power  per  unit  length  of  sheet  Is 


P ■ fw  • 
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3 2 

f »u  *r  JwB 

max 


.2wf2  -V 
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The  average  power  for  the  entire  length  of  the  Sheet  is  then 


Pt  - JP 


2 .2  2.2,  _2,2  2 .2 

■»  f Ti  Iwt  n f - B V 


IT 


6c  10''*  6c  10 
where  V is  the  volume  of  the  sheet. 
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APPENDIX  II 


SUMMARY  OF  COMPUTATIONAL  PROCEDURE  IN  CALCULATION  OF  PEAK  VOLTAGE  AND 
ADAPTATION  TO  COMPUTER  PROGRAM 

AM  relevant  permeances  are  calculated  for  the  generator  with  the  reed 
in  the  equilibrium  position.  There  are  five  such  permeances  of  which  four 
are  associated  with  the  four  air  gaps  and  the  fifth  with  the  magnetic  leakage 
paths.  The  two  gaps  at  the  fixed  end  of  the  reed  are  always  equal  and  the 
associated  permeances  are  equal  and  denoted  by  ?\  and  P2  refer  to  the 

gaps  at  the  vibrating  end  of  the  reed.  P4  is  the  permeance  of  the  composite 
of  all  leakage  paths.  P3  has  five  components  and  we  can  write  for  the  equi- 
librium condition: 

P 11  p +P  + p + p +p 

3 G3  E3  F3  L3  W3 

all  of  which  are  illustrated  in  Fig,  5 and  are  computed  by  the  following 
formulae  as  functions  of  t and  a.  All  other  dimensions  are  fixed  and 
are  the  same  as  those  of  the  sample  calculations. 


G3 


£3 


.0408 


t 


9 

.0364 


PF3  • .303 


PL3  " *2,S  ,n  V1  * 7^'- t ) * • •0916(X< 


t>]  In  t j 


and 


V*-0W‘5l"'i'*rzH) 

9 

P^  has  five  components  and  is  given  by 


V PN*  + V 4 PM4  ♦ \k  4 PW4 


PPU  m 


AH' 

V - .711 


.4  a 1.16  / 

rrr.iiv* 
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Pj  and  have  the  same  form  with  either  always  reaching  its  maximum  when 
the  other  is  at  its  minimum.  Each  has  5 components  and 


p + p + p + p + p 

Gl  El  FI  LI  KWI 


p max  p max 

Gl  . G2  . 
min  min 


.0408 


t * 2a 


El 


PE2  * .0364'' 


PFI  * PF2  * -303’ 


p max  . p max 

l!  . L2  . 
mm  min 


■2I5  V iJ 


. M70  - .0918(1  - t * 2.)3  J t 

^ A 


min  min  g 


The  total  permeance  of  the  magnetic  circuit  Pt  is  then  given  by  the  formula: 
<p, 


p m P 4 

t 4 P,  * P2  ♦ 2P, 


and  both  the  equilibrium  and  full  reed  displacement  values  of  Pt  are 
calculated.  We  use  the  equilibrium  value  P|  to  obtain  our  lo«d  line  slope 
0/H  via  the  formula 

\8 

where  /and  A are  Respectively  the  length  and  cross-sectional  area  of  the 
magnet. 

The  Intersection  of  this  load  line  with  the  demagnetization  curve  of 
Alnico  5 is  the  base  point  of  the  recoil  line  (slope  4.3  for  Alnlco  5)  along 
which  the  system  then  operates.  The  recoil  line  Is  thus  of  the  form: 

8 • 4.3  ,H  ♦ 8- 

where  is  determined  from  the  base  point.  The  maximum  value  of  ?t  corres- 
pond ing  to  full  displacement  is  then  used  in  the  relation 


8 

A H J 


•**  r max 

• - x "T 


to  obtain  the  new  load  line.  The  intersection  of  this  load  tin*  with  the 
recoil  tine  already  determined  yields  the  maximum  value  of  8 in  the  magnet. 
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The  corresponding  total  available  flux  is  then  given  by 


<*x  ■ «"“*  * 
t magnet 


and  the  flux  amplitude  in  the  reed  can  then  be  calculated  from  the  relation 
4**  « <f** 


vp.  -v 

■t  (p,  ♦ P3)(P2  ♦ p3)  ♦ %(*,  ♦ 2Pj  ♦ P2) 


The  induced  voltage  amplitude  Is  then  obtained  from 
Nm  4?*** 


.max 


10 


8 


The  computer  program  which  follows  this  procedure  is  reproduced  below 
in  Fortran  IV.  In  addition  to  calculating  the  expected  voltage  for  a given 
set  of  parameters,  the  program  also  instructs  the  computer  to  calculate 


.max 


«r 


If  the  result  is  larger  than  the  saturation  flux  density  of  the  reed  material 
(19*900  gauss),  the  calculated  voltage  is  not  attainable  and  is  marked  with 
an  asterisk  in  Table  t. 
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COMPUTER  PROGRAM  USED  TO  OBTAIN  TABLE  II 
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